The oscillating dark energy: future singularity and coincidence problem 
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lations: ideal fluid or scalar-tensor theory. It is shown that such dark energy suggests the natural 
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I. INTRODUCTION. 

It became clear recently that current universe is 
expanding with acceleration which is caused by dark 
energy (for a recent review, seeQ). However, the ob- 
servational data (for a recent review, see 01) are far 
from being complete. It is not even known what is 
the current value of the dark energy effective equa- 
tion of state (EoS) parameter w which lies close to 
— 1: it could be equal to —1 (standard ACDM cos- 
mology) , a little bit upper than — 1 (the quintessence 
dark energy Q) or less than —1 (phantom dark en- 
ergy, see [Sj for very incomplete list of related pa- 
pers). Even less is known about the evolution of 
this EoS parameter, the origin of dark energy and 
why currently the energy-density of dark energy is 
approximately equal to the one of dust matter (co- 
incidence problem). As a result, there are numerous 
proposals for dark energy which so far does not exist 
as consistent theory. 

The oscillating dark energy has been proposed in 
refs.0, IH 0- It has been shown there that such 
models very naturally resolve the coincidence prob- 
lem due to periods of acceleration and that they may 
be consistent with observations. They also serve as 
very good candidates for unification of the early time 
inflation with late time acceleration. It was argued 
that oscillated dark energy does not lead to Big Rip 



singularity. 

In the present paper the late-time cosmological 
consequences of dark energy with time-dependent, 
periodic EoS are further described. We work in 
two (mathematically equivalent) formulations: ideal 
fluid EoS description or scalar-tensor theory with 
some specific potential. The general method to 
go from one formulation to another is proposed. 
It is demonstrated for number of examples that 
such oscillating dark energy may describe the transi- 
tion from deceleration to acceleration, or from non- 
phantom to phantom era (if current dark energy is 
of phantom type). In the last case, the classical uni- 
verse may end up in the finite-time, future singular- 
ity (any of four known types of Big Rip may occur as 
is shown explicitly). The fact that oscillating dark 
energy universe may end up in the future singular- 
ity was not discussed before. It is proved that the 
oscillating dark energy may be the key for resolu- 
tion of the coincidence problem in accordance with 
earlier findings of refs. [10. Finally, we show that 
the unification of early-time inflation with late-time 
acceleration is very natural in such oscillating uni- 
verse which could be embedded in the cyclic (ex- 
panding/contracting) universe. 



II. OSCILLATING EQUATION OF STATE 
OF THE UNIVERSE. 
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Let US consider the universe filled with the ideal 
fiuid (dark energy) where the equation of state (EoS) 
parameter w depends on time t: 

p = w{t)p. (1) 
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Recent observational data (for review, see ad- 
mit the possibility of time-dependent w which lies 
currently close to —1. 

The conservation of the energy 



p + 3H{p + p) = 
and the FRW equation 



3 , 



give 



which can be integrated as 



3k2 {Jdt{l + w{t))) 
Using the Hubble rate may be found 
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H = h{t) 



3 J dt {1 + w{t)) ' 



(2) 



(3) 



(4) 



(5) 



(6) 



When w is a constant (not equal to phantom divide 
value), the standard expression is recovered 



H 



3 (1 + w) (t - t,) ■ 



(7) 



If / (1 + w{t)) = 0, H diverges, which corresponds 
to the Big Rip type singularity, which occurs when 
t = ta- In case w > — 1, for the expanding universe 
where H > 0, the cosmological time t should be 
restricted to be t > is- 

The situation where 'U'(t) is time-dependent, pe- 
riodic function may be of physical interest. Indeed, 
this can explain easily why the effective value of w 
is different at different epochs, being currently close 
to —1 (if the period of such function is comparable 
with the universe age). An illustrative example is 
given by 



w 



-1 + Wq cos LUt . 



(8) 



Here, wo,uj > 0. (Such EoS has been considered in 
ref.[5| where it was shown that it may be consistent 
with observations for some values of parameters.) 
Then Eq.® gives 



H 



2uj 



3 {wi + Wo sinujt) 



(9) 



Here wi is a constant of the integration. When 
\wi \ < wq, the denominator can vanish, which corre- 
sponds to the Big Rip singularity. When juiil > wq, 
there does not occur such a singularity. Since 



H = -- 



2uj'^Wq cos bJt 

3 (wi + Wo sinuitf' 



(10) 



when Wo cos Lut < {wo cos cut > 0) , the universe 
lives in phantom (non-phantom) phase where H > 
(ij < . In this case, the energy density p is given 

by 



3^2 (wi + wq sinwi) 



2 ' 



(11) 



which also oscillates. Hence, p > in phantom phase 
and p < in non-phantom phase. Hence, kind of 
energy transfer occurs in such oscillating dark en- 
ergy universe: the energy-density grows in phantom 
phase but decreases in non-phantom era. This may 
suggest the qualitative explanation of the fact of the 
growth of phantom energy. Moreover, in such os- 
cillating universe the unification of late-time (phan- 
tom or quintessence) acceleration with early-time, 
phantom inflation [3 naturally occurs: indeed, the 
(phantom or quintessence) acceleration epoch occurs 
periodically due to the periodic behaviour of the uni- 
verse EoS. 

Since the pressure is given by 



2H + 3H' 

from the second FRW equation, we find 
{wq COSIjjt — 1) 



3^2 (wi + Wo siiiLut) 



(12) 



(13) 



Then as long as wi + wo sin cut > 0, that is, the 
universe is expanding (H > 0), we find the following 
EoS for dark energy ideal fluid: 



2lu 



Wl 



Ky/3p 



1 + 



(14) 



Another interesting example considered in ref.| 



IS 



w = —1 + a {I + (3cosujt) 
a > , 1 > /3 > . 



(15) 



Here w > ^1 and the universe is not in phantom 
phase. Hence, this model corresponds to oscillat- 
ing dark energy where the unification of early-time 
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quintessence inflation with late-time quintessence 
acceleration occurs. Moreover, if a (1 + /3) > 1, w 
crosses w — and we find 



H 



2uj 



3a {lo (t — to) + f3 sin Lut} ' 



(16) 



Here to is a constant of the integration. In this case, 
the energy density p is given by 



P = 



iK^a^ {uj{t — to) + (3smijjt} 



(17) 



which is monotonically decreasing, oscillating func- 
tion. Since 



P 



4cj2 {{a~-l) + aPcosojt} 
in^a^ {uj{t — to) + (3 sinust} 

we find the following EoS for ideal fluid: 

p+(l- a)p 



(18) 







a(3p 



{ OLK I — 

cos( LLito -f — V 3/0 



1 



- 2(1 - a 



ap 

"2 '"^ 'i)VpY' 



2uj ^ 
2a ~ 1) p2 



(19) 



In a more realistic situation, universe is filled with 
matter and dark energy. The matter, in general, 
interacts with the dark energy. In such a case, the 
total energy density ptot consists of the contributions 
from the dark energy and the matter: ptot = p + Pm- 
If we define, however, the matter energy density p^ 
properly, we can also define the matter pressure pm 
and the dark energy pressure p by 



Pn 



-Pn 



P 

3H 



P = Ptot - Pn 



(20) 



Here ptot is the total pressure. Hence, the matter 
and dark energy satisfy the energy conservation laws 
separately, 

p„j + 3H (p„i + p,„) = , p + 3H {p + p) ^ . 

(21) 

We now consider the case that universe is filled 
with the oscillating dark energy and the matter ideal 
fluid with a constant EoS parameter Wm- Pm = 
WmPm- Since Wm is a constant, one finds 

Pm = Poa-3(^+'"'") . (22) 

Here po is a constant. FRW equation leads to 



P = -^H - poa 



(23) 



Using the conservation law |2Jl, we find 



-i{l + w{t)) = -\np 

_ -^HH + 3(1 + w„,) Hpoa-^^^+^'^'> 



(24) 



That is, 



Note that the denominator of the r.h.s. in H25|l is 
always positive since the denominator is nothing but 
energy-density. 

Let us consider the transition from non-phantom 
to the phantom era, where H = 0. Then 1 + u'(i) < 
as long as Wm > —1. Thus, in order that the 
transition occurs, w(t) should become less than — 1 
before the transition. 

As an example, we consider the scale factor 



a{t) = ao 



t 



t,-t 



(26) 



with constants ao, tg, and ho > 0. Ea. (|26|l gives, 
„ hots 

— — 7 7 5 

t {ts - t) 

a _ hpts {2t - (1 - ho) ts) 
a " t^ {ts - t)^ ■ 



(27) 



The transition from non-phantom to phantom era 
occurs at t — ts/2, and if ho < 1, the transition 
from decceleration to acceleration epoch occurs at 
t = Ml^iM ^ Vsing (|23l, w{t) corresponding to ^ 
is given by 



w{t) = -1 

_ [ 2Ms i2t~ts) 
\ nH^ (ts - t)^ 

+ {l + Wm)poao 



ts-t 



^2*2 (ts - ty 



-Poa-o 



3(l+to„) 



ts~t 



-3(l+iD„)/io ■ 



(28) 
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Now the energy-density is given by 
Chit's 



Pit) 



3(1 + M)„) 



t 



- 1 



(29) 



Since the second term in is negative, p{t) is not 
always positive. The negative energy density is, of 
course, not physical. When t ^ 0, if 3/io (1 + Wm) > 
2, the second term dominates and the energy density 
becomes negative. On the other hand, when t — > 
ts, if 3/io (1 + Wm) < —2, one finds that p becomes 
negative, again. Hence, physically allowed region of 
the parameters could be 



-2 < 3ha (1 + w„0 < 2 



(30) 



For simplicity, we consider the case SHq (1 + Wm) 
2. Then the expression H29|l reduces to 



Pit) = ^QW , 



Qit) 



01,2^2 2 

"2 - poao its ~ t) ■ (31) 



(t - t,Y 



Since Q(t) is monotonicaly increasing function as 
long as i < ts, if 



g(o) 



3/1^ 



(32) 



p(t) is always positive. 

From the second FRW equation, the pressure may 
be found: 



1 

P = 2 



2R + 3H^] - WmPoa 



Hence, 



Pit) 



2hots {2t-ts)+3hlt 



-3(l+u)„0 



2+2 



(33) 



kH^ {ts - t)" 



-WmPoaQ 



3(l+-!«„) 



ts-t 



-3h„{l+w,„) 



.(34) 



At the phantom-non-phantom transition point t = 
ts/2, one gets 

Pit) + Pit) 

= - (1 + Wm) Poa^ 



- 1 



(35) 



which is negative, corresponding to w{t) < —1 there. 

As the next example, the following scale factor 
may be considered: 



a = aoe'"'*+^'''"'^* 



(36) 



with constants ho > hi > and w > 0. Since 

H — ho + hi cosLut . (37) 

Here, the universe lives in non-phantom 
(phantom) era when 27rn < ut < 2Tm + tt 
{2Tm — TT < ujt < 27to), where n is an integer. The 
corresponding, time-dependent EoS parameter w(i) 
is found to be 

w{t) = -1 

( 2hi . 

— < — smuji 

(1 + Wm) poao-'('+'"")e-3(i+-'")(''o*+^ 

f 3 2 

X < — (ho + hi cos ujt) 

Now the energy density p and pressure p are 

3 2 
pit) = ^ (/iQ + hi cos ujt) 

p{t) = — |2/iiw sinuit — 3 {ho + hi cosa;t)^| 



-WmPoGo 



-3(l+tii,„) _3(l+to„)(/iot+iLL sin Lut) 



(39) 



The p{t) H39|l becomes negative for large negative t. 
If we restrict the parameters to satisfy 

^ {ho - hi)^ - poao > , (40) 

it follows p{t) > when t > 0. 

It is interesting to understand if the oscillating 
dark energy may bring the universe evolution to Big 
Rip singularity. (One should not forget that Big 
Rip is typically classical effect which seems to dis- 
appear with proper account of quantum corrections 
IE El-) The classification of future, finite-time 
singularities is given as|ld|: 

• Type I ("Big Rip") : For t i^, a ^ cx), 
p —^ 00 and \p\ —^ 00 

• Type II ("sudden") : For t ^ t^, a as, 

p ^ Ps or and \p\ 00 
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Type III : For t ^ ts, a 

\p\ -> oo 



oo and when a 7^ 1 and for a = 1 



• Type IV : For t — > t^, a — > Os, p — s- 0, |p| — > 
and higher derivatives of H diverge. This 
also includes the case when p (p) or both of 
them tend to some finite values while higher 
derivatives of H diverge. 

Here ts, and ps are constants with ^ 0. 
Let assume 

H^ho ~ ty" , (ho > 0) , (41) 

near the singularity t ^ t^. Then one gets 



H 







{a = 0) 



Since p cx \/H and p oc iH^ + 2H, in terms of a, 
the above four types of the singularities could be 
classified as 

Type I: a > 1, Type II: > a > -1, Type III: 
1 > a > 0, Type IV: a < — 1 but a is not an inte- 
ger (for account of Hubble rate dependent terms in 
inhomogeneous EoS, see [lT|'). 

Without the matter, in terms of w{t), by compar- 
ing (gll) with iP, we find 



l + w{t) - {t.-tf 



(43) 



Thus, all four types of future singularities are real- 
ized here for above values of parameter a. 

It is interesting to study what happens with the 
classical future singularity in the presence of matter. 
For this purpose, we assume Ea. (|41|l . Then from 
(EH, we find 

l + w{t) 



f '2a , , 
^/io (ts-t) 



(44) 



1 -I- w{t) 



+3il + w„^)poao (i, - t)-''°(i+'"'")) 

\ Hi 

-poaoiU-t)-''"^'+''-A'\ (45) 



In case a > 1, which corresponds to Type I sin- 
gularity, one obtains 

l + w{t)^^iU-tr-'^0. (46) 

if 1 + w„i < 0. The case 1 -I- Wm > is excluded 
since the denominator of the r.h.s. in H25|) becomes 
negative. 

In case a = 1, which also corresponds to Type I 
singularity, 

{-2§/lO+3(l + -l«„OPQao -r , /-, , N r, 

— j-Ta it /in 1 -f Wjj, — 2 

if /lO (1 + w™) = 2 
(47) 

The case ho {1 + Wm) > 2 is also prohibited since the 
denominator of the r.h.s. in H25|) becomes negative. 

In case < a < 1, which corresponds to Type HI 
singularity, we find 



2a , 

1 + ^(i) ^—{ts- t) 



\a-l 



3ho 



00 , 



and in case a = 0, 



l + w{t) 



3(1 



I POCIQ 



r/lg — PqUq 



(48) 



(49) 



The case —1 < a < 0, which corresponds to 
Type II singularity, and the case a < —1, which 
corresponds to Type IV, are excluded since the de- 
nominator of the r.h.s. in H25II becomes negative. 
It is interesting that the Type II and IV singulari- 
ties are excluded with the account of matter, while 
they occur without matter. Thus, we demonstrated 
that within oscillating dark energy universe, all four 
types of future singularities may occur on the clas- 
sical level. This is our main qualitative result: the 
oscillating dark energy of general form does not pre- 
vent the universe to end up in the future singularity 
of any from four known types. 
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III. SCALAR-TENSOR DESCRIPTION AND 
COINCIDENCE PROBLEM 

In this section we will present the oscillating 
dark energy with time-dependent, explicit EoS in 
an equivalent, scalar-tensor description following to 
the method developed in ref.l2]. (Note that such 
method may be applied to include also the time- 
dependent bulk viscosity |l3l introduced for dark 
energy description in ref.jlJI-) Let us start with the 
following action 



S = 



-v{<t>)} , 



(50) 



with an adequate function f{4>). The following so- 
lution of FRW equation exists (151]: 

= t, H = f{t). (51) 

Then by comparing (|51|l with (jSJ , one may identify 

fit) = hit) . (52) 

For the action H5()|l . the explicit EoS is given by: 



Comparing H41|) with (|52|l . one also finds that in 
scalar-tensor theory (|Sn|l . the four types of future 
singularities could be realized by the choice 

/(0)~/o(t.-0)-" = /o(/^-" , (56) 

with positive constant /q. Here Lp = tg — cj). In this 
case, we obtain 



^i^) - its 



(57) 



^(<^) = ;^(3/o'(is-0)~'"+«/o(t.-0)-"-') , 

One should bear in mind that even FRW equations 
are the same in EoS or scalar-tensor descrition, the 
emerging universes may be different. Indeed, the 
number and stability of FRW solutions may not co- 
incide, the newton law may be slightly different, etc. 

As an extension of 1)5 Of) , the matter with constant 
EoS parameter w — Wm may be added into the 
action ((50|l . By using a single function git), if we 
choose the scalar potentials in the action (|50|) with 
matter as 



n^) = -(3g'(0)2-fg"(0)) 



UJ„, ~ 1 



-5oe 



(58) 



4/' 



p 



(53) 



the following solution of FRW equation may be con- 
structed (see also 0): 



Thus, we demonstrated how an arbitrary oscilat- 
ing dark energy may be rewritten explicitly in 
mathematically-equivalent scalar-tensor form. 

Hence, for Eq.®, the following equivalent de- 
scription in scalar-tensor form may be constructed: 



ni<t>) 
vi<k) 



'iuj'^wo COS uj4> 
(wi -I- Wo sinwt/ 
(4 — 2wo cosLuq 
(wi -|- wq sinujq 



(54) 



For EoS l|15ll . the equivalent scalar-tensor theory de- 
scription 12] leads to 



Viqb) 



4u;^ (1 + /3 cosuff 



3aK^ {^i(f> — to) + Psinujcf)} 
[2 - a - a(i cos ^0) 
'6d^K? {a;((^ — to) + sinw^}^ 



(55) 



t , g'it) , 



a = aoe^(*) , ao ^ \^ ) ) (59) 



In the case of scale factor (|26|) 



;(0) - In 



tf. 



ho 



and in the case of 



gi(j)) = In ( ho(p H — - sincj^ 



(60) 



(61) 



Let us now consider a little bit different function 



gi(f>) — a\n-!^ got + — cos uj(f> 



(62) 
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Here a, go, Qi, and co are positive constants. Then 
the Hubble rate is given by 



H 



« (go - 91 sinut) 
got + — cos 



(63) 



li t > 1/u!, the Hubble rate H is always positive and 
if go is large enough compared with gi and uj, we 
find H < 0. Taking p and p in Q as total energy 
density and pressure, w{t) is given by 



w{t) 



2 

+ — 



251WCOS ut 



3ck 3a {go — §1 sin u!t)^ 
When go ^ gi, the expression H64I) reduces to 

2giti;cos tot 



wit) 



^ 2 
3a 



Then if 



> II 



(64) 



(65) 



(66) 



w oscillates and crosses w = —1/3. Hence, the uni- 
verse expiriences the transition from the deccelera- 
tion to the acceleration epoch. 

Having in mind the oscillating nature of dark en- 
ergy under consideration, it is expected ^lEIll that 
such type of dark energy may naturally solve the co- 
incidence problem (for recent discussion and list of 
refs., see |l6j). Let us confirm this in both formula- 
tions of the oscillating dark energy under consider- 
ation. We now separate ptot into the contributions 
from the dark energy and matter: ptot = P + Pm- 
If the EoS parameter Wm of the matter is constant 
as in 1)22(1 . by using the first FRW equation one 
gets 

r = — = -1 H ^ . (67) 

Pm ^ PrnO 



For Eq.^, it follows 



r = -1 



3a a] 



2 3{l+w^) 







(50 - 51 smujt) 



s 3(l+«;,„)a-2 

X I ant -I cos utj . (68) 



With the choice a = 2/3{l + w,n), the Eq.lgHIl re- 
duces to 



r = -1 



3a a, 



2„3(1-I-M„) 



k'^PtuO 



[go - gi sinujt) . (69) 



Then r has a maximum and minimum as 



3/ 
a a, 



2 3{1+Wra) 



< -1 



3a a. 



2„3(l-l-ii 







k'^PtuO 



[go - gi) < r 



(30 + 51) 



(70) 



Hence, if -1 + {3a^af^^+'"^^ '-^ 



/t^Pmo) (50 ±51) ~ 

0(1), the ratio between the energy densities of 
the dark energy and matter could be always 0(1). 
Hence, oscillating dark energy may solve the coinci- 
dence problem. 

In case of or from one gets 



r = -1 



K2pot2-3ho(l+™™) 



(71) 

We now assume (ISUt . 

When t ts, if 3ho {1 + Wm) ^ —2, r diverges 
but if 3ho (1 -|- Wm) — —2, r is finite: 



-1 



(72) 



On the other hand, when < = 0, if 3ho (1 -I- w„i) 7^ 2, 
r diverges but if 3ho (1 + Wm) = 2, r is finite: 

r -> ro . (73) 

If we write t = jts, when 3ho (1 -I- Wm) — 2, wc find 



n 



n^pot^^ (1 - 7)' 



(74) 



In the present universe, < 7 < 1 and 7 '--^ 0(1). In 

order that tq - 0(1), 3/igao'^+"''"V'«^Pois should be 
larger than 1 and 0(1). Then we also find r; ~ 0(1). 
This tells that if the ratio of the energy densities of 
the matter and the dark energy is 0(1), even in the 
present universe, the ratio is 0(1), which may solve 
the coincidence problem. 

For the model in (|36() or H61() . by using H39|l one 
obtains 



r = -1 



3a, 







K^PO 



(ho + hi cos wt) 



xe 



3{l+Wm){hot+!^ sinwt) 



(75) 



which increases exponentially with time if Wm > — 1 . 

Thus, we confirmed that oscillating dark energy 
may naturally solve the coincidence problem in both 
(EoS or scalar-tensor) formulations which is not 
clear from the beginning because the mathematical 
equivalence does not mean the physical equivalence. 
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IV. DISCUSSION 

In summary, we presented two equivalent formu- 
lations of oscillating dark energy. It is demon- 
strated that such dark energy may describe very 
naturally most relevant late-time cosmological phe- 
nomena: transition from deceleration to acceleration 
era and transition from non-phantom to phantom 
epoch (if the universe currently enters such epoch). 
If the universe currently lives in phantom era, it is 
demonstrated explicitly how the oscillating dark en- 
ergy may bring the evolving universe to the one of 
four known types of finite-time, future singularity 
(on classical level). Moreover, oscillating dark en- 
ergy may suggest nice resolution of the coincidence 
problem as was suggested in refs. HH. Finally, such 
dark energy unifies in a clear and natural way the 
early-time inflation and late-time acceleration what- 
ever (quintessence, cosmological constant or phan- 
tom) nature they have. However, in order to check 
if such unification is realistic one, it is necessary to 
analyse the details of inflation, especially, the pre- 
heating scenario. This lies beoynd the scopes of our 
paper because some work in this direction (for phan- 
tom inflation) was already done in 0. It is also in- 
teresting that using method [l3 the oscillating dark 
energy may be presented as some special form of 
modified gravity (for recent review, see 01 )• 

One can go further and propose the following cos- 
mological model: 

H = ho sin ujt {I + hi sin rriLut) . (76) 

Here m is a positive integer, hg, hi, to are pos- 
itive parameters and it is assumed that m 3> 1, 
hi < 1. When 2?™ < utt < 27m + tt (n:integer), 
H is positive and universe is expanding but when 
2Tm — TT < Lut < 2Tm (n:integer), H is negative and 
universe is contracting. Hence, the universe is os- 



cillating with period 2-k/u!. The conjecture could 
be that our cyclic universe (compare with jlSi] ) lives 
currently in the expanding era (repeating the accel- 
eration/superacceleration circles each 14-15 billion 
years) while the transition from accelerating to con- 
tracting era occurs with much bigger period of time. 
Since it is assumed m 3> 1 , adiabatically one obtains 

H = hauj {cos ojt (\ + hi sin mujt) 
-\-himsin ut cos mcut} 
= mwhohisinujtcosmujt . (77) 

Then in the expanding phase (sin a; > 0), if 2?™ — 
7r/2 < w < 2?™ + 7t/2 (n:integer), we find H > Q, 
that is, the universe enjoys the phantom era and if 
2nn + 77/2 < Lo < 27rn -I- 3n/2 (n:integer), H < 0, 
that is, the universe is in non-phantom era. In the 
contracting phase (sinw < 0), if 27rn — 7r/2 < w < 
27rn-t-7r/2 (n:integer), H < Q and if 27Tn+Tr/2 < uj < 
2?™ -I- 37r/2 (n:integer), H > 0. Hence, the universe 
iterates phantom and non-phantom eras with the pe- 
riod 217 jniLo. There is no problem to construct the 
EoS and/or scalar-tensor description corresponding 
to scale factor (|76|l . Nevertheless, the existing data 
cannot confirm/reject cyclic cosmological models at 
high confidence level so one awaits the next gener- 
ation observational data which will be able to re- 
construct the realistic evolution of the cosmological 
parameters. 
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